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Exoplanets at small orbital distances from their host stars are submitted to intense levels of energetic radiations, X-rays and extreme 
ultraviolet (EUV). Depending on the masses and densities of the planets and on the atmospheric heating efficiencies, the stellar 
energetic inputs can lead to atmospheric mass loss. These evaporation processes are observable in the ultraviolet during planetary 
transits. The aim of the present work is to quantify the mass-loss rates (in), heating efficiencies (if), and lifetimes for the whole 
sample of transiting exoplanets, now including hot jupiters, hot neptunes, and hot super-earths. The mass-loss rates and lifetimes are 
estimated from an "energy diagram" for exoplanets, which compares the planet gravitational potential energy to the stellar X/EUV 
energy deposited in the atmosphere. We estimate the mass-loss rates of all detected transiting planets to be within 10 6 to 10 13 g s" 1 
for various conservative assumptions. High heating efficiencies would imply that hot exoplanets such the gas giants WASP-12b and 
WASP-17b could be completely evaporated within 1 Gyr. We further show that the heating efficiency can be constrained when in is 
inferred from observations and the stellar X/EUV luminosity is known. This leads us to suggest that r] ~ 100% in the atmosphere 
of the hot jupiter HD 209458b, while it could be lower for HD 189733b. Simultaneous observations of transits in the ultraviolet and 
X-rays are necessary to further constrain the exospheric properties of exoplanets. 
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1. Introduction 

More than 500 extrasolar planets have been detected so far. A 
large fraction (~ 30%) of them lie at close distances from their 
host stars, below 0.1 au. Early on, the fact that most detected 
planets were so close to their stars and massive enough to be 
most certainly gaseous, deeply questionned theoreticians about 
their origins and their fates. "Hot jupiters" cannot form at their 
observed locations, but have to migrate inward (Lin et al. 1996). 
Meanwhile, why did these planets stop migrating at such dis- 
tances is still an open question. In contrast, one could wonder 
what happened to those giant planets that migrated further in. 
Of particular interest, is the question of the atmospheric stabil- 
ity at extreme levels of irradiation. How fast can giant planets 
be despoiled of their atmospheres? Couldn't lighter, Neptune- 
or Earth-mass planets form or migrate as well into environ- 
ments that would questioned their atmospheric stability? Can 
such "evaporation" of exoplanets explain the apparent desert of 
planets below ~ 0.01 au? In fact, the closest planets should also 
be the easiest to detect because the transit frequency and proba- 
bility increase inversely to the semi-major axis. 

Observations of exoplanet transits hold the keys to these 
questions. For the > 100 known transiting planets, it is possible 
to measure the planet-to-star radius ratio and, for the few planets 
observed in the ultraviolet, constrain the size and mass-loss rate 
of their evaporating upper atmospheres. Transit observations of 
hot jupiters in the stellar Lyman-a' emission of neutral hydrogen 
(Hi Lya at 1 215.67 A) yield typical estimations of the mass- 
loss rates of m ~ 10 10 to 10 11 g s _1 for the planets HD 209458b 
(Vidal-Madjar et al. 2003, 2004, 2008; Ehrenreich et al. 2008; 
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Linsky et al. 2010), HD 189733b (Lecavelier des Etangs et al. 
2010), and WASP-12b (Fossati et al. 2010). Such values are pre- 
dicted by numerous theoretical works and imply that the evo- 
lution of known hot jupiters is not significantly impacted by 
atmospheric evaporation (Lammer et al. 2003; Lecavelier des 
Etangs et al. 2004; Baraffe et al. 2004, 2005, 2006; Yelle 2004, 
2006; Jaritz et al. 2005; Tian et al. 2005; Garcfa-Munoz 2007; 
Holmstrom et al. 2008; Stone & Proga 2009; Murray-Clay et al. 
2009). 

Lecavelier des Etangs (2007) compares the stellar energy re- 
ceived by the upper atmospheres to the gravitational energies 
of the planets in a so-called "energy diagram". This diagram 
allows estimations of the mass-loss rates and lifetimes for the 
population of close-in planets. This author uses the sample of 
exoplanets known as of June 2006, including ten transiting plan- 
ets. Recently, Davis & Wheatley (2009) have provided an up- 
dated energy diagram containing 36 transiting hot jupiters, while 
Lammer et al. (2009) have determined the mass loss limit for 57 
transiting planets. These studies pointed out that there should ex- 
ist an evaporation-forbidden region into which a lost population 
of planets have been significantly eroded by this process, pos- 
sibly leaving evaporation remnants, coined "chthonian planets" 
(Lecavelier des Etangs et al. 2004). 

In this article, we aim at providing comparative estimations 
of the mass-loss rates and lifetimes of the whole sample of tran- 
siting planets currently known in October 2010 (105 planets). 
As the sample of transiting planets steadily increases in time, 
the explored ranges of mass and density expand to new kinds of 
exoplanets: hot neptunes and hot super-earths. The existence of 
such close-in low-mass and high-density planets is a test for at- 
mospheric evaporation theories. Moreover, predictions from the 
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energy diagram enable to estimate the observable transit signa- 
ture of evaporating planets (e.g., Ehrenreich et al. 201 1). 

The method to derive the energy diagram is detailed in 
Sect. |2] Afore-mentioned theoretical studies point out the role of 
X and extreme ultraviolet (EUV) stellar radiations as the main 
source of exospheric heating leading to mass loss. As recent 
X/EUV surveys were dedicated to exoplanet host stars (Kashyap 
et al. 2008; Sanz-Forcada et al. 2010; Poppenhaeger et al. 2010), 
we examine the impact of the constraints brought by these stud- 
ies on the energy diagram (Sect. |3). In a few cases, measure- 
ments of the stellar X flux and estimations of the mass-loss rate 
from observations of ultraviolet transits make it possible to con- 
strain the mechanisms of atmospheric evaporation in terms of 
heating efficiency. This is discussed in Sect. [4] Finally, the esti- 
mation of the atmospheric mass-loss rates raises questions about 
the atmospheric stability and the lifetime of evaporating planets, 
which is discussed in Sect. [5] 

2. Updated energy diagram 

The energy diagram for extrasolar planets measures the gravi- 
tational potential energy per mass unit of a planet, as a func- 
tion of the X and EUV irradiations received by the planet. Next, 
we recall the main steps of the calculations further detailed in 
Lecavelier des Etangs (2007). All the properties of the plan- 
ets and host stars used for these calculations are extracted from 
the Extrasolar Planets Encyclopaedia (Schneider 2010) and are 
given in Table Q] 

The framework is that of an energy-limited escape, i.e., an at- 
mospheric mass unit can escape when the gravitational potential 
(dE' p ) is filled by the stellar X/EUV energetic input (d£x/Euv) 
into the planet upper atmosphere, 

- dE' p = 77d£x/Euv- (1) 

The key factor r\ is the heating efficiency in the exoplanet 
thermosphere that encapsulates most of the physics of the prob- 
lem. This factor value is high when most of the incoming stellar 
energy is effectively used to escape the atmosphere. It depends 
mainly on the fraction of energetic photons absorbed by atmo- 
spheric atoms, how deep in the atmosphere the absorption takes 
place, and how thick the absorbing region is. 

Lecavelier des Etangs (2007) consider the extreme case 
r\ - 1, where all the stellar flux is used to escape the at- 
mosphere. Considering (without quantifying) energetic losses 
due to thermal emission by atmospheric hydrogen, Tian et al. 
(2005) model the atmospheric escape process of the hot jupiter 
HD 209458b using r\ - 0.15. This value was initially chosen 
by Watson, Donahue & Walker (1981) in their pioneering work 
about Earth's atmospheric escape. A similar value (77 = 0.1) 
was also employed by Valencia et al. (2010), who estimate the 
atmospheric loss for Corot-7b. In fact, the value of 77 for the 
strongly irradiated atmospheres of extrasolar planets remains 
largely unconstrained (though see Cecchi-Pestellini et al. 2009 
and Sect.|4j>. In the following, we treat 77 as a free parameter with 
a possible value between 0.01 and 1, until we detail in Sect. [4] 
how it could be possible to obtain observing constraints on this 
parameter. 

We present the updated energy diagram in Fig. [1] with 98 
transiting planets. The vertical axis of the energy diagram mea- 
sures the potential energy required to escape the planetary atmo- 
sphere. The potential energy per atmospheric mass unit dE' p /dm 
of a planet of radius R p and mass M p includes the contribution of 
tidal forces and is calculated according to Eq. (8) in Lecavelier 
des Etangs (2007). 
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Fig. 1. The energy diagram for 98 transiting planets: the energy 
needed to escape a unit of mass of the planet atmosphere ver- 
sus the X/EUV flux reaching the top of the atmosphere per unit 
of time. The dotted lines indicates constant mass loss rates of 
m = 10 14 to 10 5 g s -1 (from left to right). Data points are calcu- 
lated with 77 = 15%. Variations of 77 between 1% and 100% are 
represented with horizontal gray error bars. Coloured points and 
labels indicate transiting planets of particular interest which are 
discussed in the text. 



The horizontal axis of the energy diagram measures the 
amount of stellar energy available on top of the planetary atmo- 
sphere. Cecchi-Pestellini et al. (2009) have studied how EUV 
photons (30-912 A) and X-rays (1-30 A) can heat the upper and 
lower atmospheric layers, respectively. This heating makes the 
atoms and ions of the upper atmosphere escaping the gravita- 
tional well of the planet. Lecavelier des Etangs (2007) estimates 
the amount of X/EUV energy received per unit of time as 



d£x/Euv _ 1 2 -2 j 

J t - 1 R P a p L X/EUV- 



(2) 



The X/EUV luminosity Lx/euv (1-912 A) is generally not 
known for host stars of transiting planets (but see Sect. |3). 
Lecavelier des Etangs (2007) uses the Wood et al. (1994) corre- 
lation between the stellar rotation velocity v lot and the EUV flux 
^euv measured in the S2 bandpass of the Rosat satellite (110- 
195 A) as a proxy to estimate the EUV flux of the 11 transited 
stars known in 2007, 



F EU y(l au) = 4.6 



Vrot 



2.0 km s 



1.4 



erg cm 2 s 1 . 



(3) 



The luminosity and flux at 1 au are related through Leuv = 
47r(l au) 2 Feuv(1 au). We caution here that such an activity- 
rotation relation is only valid for stars with a convective zone. 
We therefore consider this estimation to be doubtful for stars 
earlier than ~F5 (about 1 1% of transited stars). In particular, the 
position of WASP-33b (A5 star) in Fig. [TJ should be regarded 
with caution. 

The measured v, ot sin 7* have been retrieved from the liter- 
ature for 98 stars in Table [TJ We also assumed that the sine of 
the stellar inclination 7* is ~ 1 for transiting systems. The data 
points in Fig. [TJ are calculated assuming 7/ = 0.15. The value 
of this factor dominates the uncertainties on such calculations. 
Therefore, the positions of planets in the diagram are shown 
with conservative error bars representing possible 77 values in the 
range [0.01, 1]. Curves of constant mass-loss rates m are over- 
plotted on Fig. [TJ All transiting planets in the diagram have m 
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Fig. 2. Mass-loss rate (77 = 0. 15) as a function of the planet mass 
for 98 transiting planets. The vertical error bars represent a vari- 
ation of 77 from 0.01 (lower values) to 1 (higher values). Planets 
of interest are coloured accordingly to Fig. [T] 

values between 10 6 g s _1 and 10 13 g s _1 , assuming all possible 
values of 77. 

Two 'outliers', WASP-33b and HD 80606b, appear in the di- 
agram with the highest and lowest mass-loss rates, respectively. 
Meanwhile, the mass-loss rate of WASP-33b is probably over- 
estimated for the afore-mentioned reason. The mass-loss rate of 
the extremely eccentric (<? = 0.93) planet HD 80606b is likely 
underestimated, at least when the planet is at periastron, since 
the mass loss rates are calculated with the semi-major axis of 
transiting planets instead of the periastron distances. 

The iso-mass-loss rates are represented by linear curves in 
Fig. Q] The mass-loss rate is also plotted as a function of the 
planet's mass in Fig. [2] Both Figs. Q]and[2] show that maximum 
mass loss rates (10.5 $ logm $ 12) are reached with 77 = 0.15 
for planets with masses between ~ 0.4 and ~ 4 M^. In Fig. [21 
there is a dearth of low-mass planets with high m, which is ex- 
pected if low-mass planets are despoiled of their atmosphere at 
high m. The remaining cores (or chthonian planets) would be 
more difficult to detect in transit. Assuming a different value for 
77 would simply scale up or down the resulting mass loss rates, if 
one assumes a similar value of 77 for all planets. See Sect.|4]for a 
discussion about this aspect. 

3. Measurements of X/EUV fluxes 

To take into account the highest-energy photons able to heat the 
upper atmospheres of the transiting planets, it is necessary to es- 
timate the X/EUV fluxes of planet-hosting stars. Observational 
and archival surveys have been started for this purpose. Kashyap 
et al. (2008) carried out a survey of X-ray emission combining 
archival and targeted data from Asca, Exosat, Einstein, Rosat, 
XMM-Newton, and Chandra missions. These authors present 
X-ray flux values or upper limits for 235 planet-hosting stars. 
These values cover a pass band of 0. 1—4.5 keV (3-124 A). There 
are 35 transiting systems among their sample, among which 
only 5 have a detected X-ray emission (GJ 436, HD 209458, 
HD 189733, OGLE-TR-10, and SWEEPS-11). 

Sanz-Forcada et al. (2010) study a sample of 65 planetary 
systems with XMM-Newton and Chandra archival data. The data 
sets are the same as in Kashyap et al. (2008) but some estima- 
tions of stellar X luminosities are revised. According to Sanz- 



Fig. 3. The energy diagram for 35 transiting planets with X-ray 
emission detections or upper limits (arrows) reported in the lit- 
erature. Caption is the same as for Fig. Q] Data points are calcu- 
lated with 77 = 15%. Planets of interest are outlined with different 
colours. 



Forcada et al. (2010), these differences appear in some cases 
(noticeably for HD 209458) where the stellar proper motion is 
not taken into account, leading to suspicious source identifica- 
tion. Concerning transiting systems, Sanz-Forcada et al. (2010) 
report lower values than Kashyap et al. (2008) for HD 189733 
(by 0.2 dex) and GJ 436 (by 1 dex). They give an upper limit 
for HD 209458, significantly lower (by 0.7 dex) than the values 
reported elsewhere. 

Poppenhaeger et al. (2010) present a similar work as 
the above-quoted studies. They add new XMM-Newton data 
points for stars previously devoided of X-ray measurements. 
Unfortunately, their new data set does not encompass additional 
transiting systems. They report similar X-ray luminosity values 
for HD 189733 and GJ 436, although their calculated uncertain- 
ties are larger than those estimated in Kashyap et al. (2008). 

For some active stars, the X-ray luminosity may vary outside 
the range allowed by the uncertainties reported in these works. In 
fact, a detailed analysis of archival XMM-Newton observations 
of the active star HD 189733 at two different epochs show that 
these uncertainties probably underestimate the intrinsic source 
variability (Pillitteri et al. 2010). Pillitteri et al. (2010) indeed 
find that the X-ray flux of this star has variated by 45% between 
2007 and 2009. 

The values of Lx/euv found in the literature have been in- 
jected in Eq. (f2]i to obtain the energy diagram shown in Fig. [3] 
Among the five transiting systems with measured X-ray emis- 
sions (Kashyap et al. 2008), HD 189733b and GJ 436b stay at 
about the same location as in Fig. Q]for 77 = 0.15, with m be- 
tween 10 10 and 10 n g s 1 for HD 189733b and ~ 10 9 g s 1 
for GJ 436b. In contrast, the mass loss rate of HD 209458b is 
about 1.5 order-of-magnitude lower than with the previous es- 
timation. Noticeably, the m value for OGLE-TR-lOb dramati- 
cally increases by ~ 2 orders of magnitude. Finally, the case of 
SWEEPS-11 is troublesome. There is no reported value for the 
rotational velocity of this star in the literature, so we were not 
able to include this object in Fig.Q] The X/EUV flux reported by 
Kashyap et al. (2008) for this star is more than 10 times higher 
than for HD 189733. SWEEPS-11 is a 20 th -magnitude star be- 
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Table 2. Multiple X-ray detection reports in the literature for stars hosting transiting planets 



Star 




log(L x [ergs ']) 








Hiinsch et al. 


Kashyap et al. 


Sanz-Forcada et al. 


Poppenhaeger et al. 


Pillitteri et al. 




(1999) 


(2008) 


(2010) 


(2010) 


(2010) 


HD 189733 


28.44" 


28.43 ± 0.01" 


28.18" 


28.26 ±0.12" 


27.72 6 -28.29 6 


HD 209458 




27.02 ± 0.20* 


< 26. 12* 






GJ 436 


26.85" 


27.16 + 0.19" 


25.96 6 


27.16 ±0.34" 





Notes. " Rosat/PSPC. h XMM-NewtonfEPlC. 



longing to a crowded region in the Galactic bulgeQ The possi- 
bility that the detected emission does not come from this star is 
therefore not negligible. 

Overall, there are multiple reports of the X-ray emissions 
of only three transiting systems in the literature. These are 
HD 189733, HD 209458, and GJ 436. The different values re- 
ported for these three systems are summarized in Table 

4. Constraints on the heating efficiency 

For the majority of transiting planets, the two largest uncertain- 
ties in the calculation of m are the heating efficiency 77 and the 
X/EUV luminosity. Cecchi-Pestellini et al. (2009) have stud- 
ied the relative role of X/EUV radiation in the heating of the 
hydrogen-rich atmosphere of HD 209458b. In their model, the 
flux of stellar X/EUV photons incident upon the planetary atmo- 
sphere of solar-like composition photo-ionizes the gas produc- 
ing a flux of high energy photo-electrons, which deposit their 
energy into the gas. They have calculated heating efficiencies as 
a function of photon energy, electron fraction in the gas, and ver- 
tical column density of the atmosphere. The results strongly de- 
pend on the coupled ionization and density structures as shown 
by Yelle (2004), Garcfa-Munoz (2007), and Murray-Clay et al. 
(2009). For a fractional electron concentration of 10%, Cecchi- 
Pestellini et al. (2009) found heating efficiencies of ~ 0.9, 0.75, 
and 0.7 for column densities 10 19 , 10 21 , and 10 22 crrT 2 , re- 
spectively, and for photon energies of 50 eV (A = 248 A), 300 e V 
(41 A), and 1 000 eV (12 A), respectively. These efficiencies de- 
crease to 0.18, 0.04, and 0.02 when the electron fraction is 10~ 6 . 

The work of Cecchi-Pestellini et al. (2009) shows that de- 
pending on the atmospheric properties, the heating efficiency can 
take nearly all possible values. Reversly, measuring the heating 
efficiency 77 could provide some constraints on these atmospheric 
properties. In the few cases where transits have been observed in 
the UV, it is possible to constrain m. This can be done (i) numer- 
ically using a particle simulation with m set as a free parameter 
to reproduce the absorption during transit caused by the cloud of 
escaping hydrogen atoms (Lecavelier des Etangs et al. 2010), or 
(ii) analytically by assuming a density profile - usually a power- 
law - for the hydrodynamically escaping atmosphere (Linsky et 
al. 2010). If the X/EUV luminosity is known, UV transits can 
thus bring constraints on 77. In fact, starting from Eq. ([!}, the 
heating efficiency can be expressed as the ratio of the two pow- 
ers Lx/euv and L hl depending on the star and planet properties, 
respectively, 

I = 7 ■ (4) 

i-X/EUV 

If tidal forces are neglected, the "mass loss power" L„ has a 
simple expression, 

L ri1 = (l6n/3)Gmp p a 2 p , (5) 
1 This comment could also apply to OGLE-TR-10. 



where p p = M p /(4/3jrR p ) is the mean density of the planet. 

In the following, all the calculations are performed taking 
tidal forces into account. To constrain the value of 77, we intro- 
duce in Fig. |4] a "power diagram" for extrasolar planets. Only 
two transiting planets have both observational constraints on m 
and Lx/euv: HD 209458b and HD 189733b. For HD 209458b, 
we have considered the value of 10 10 g s for the mass-loss 
rate given by Vidal-Madjar et al. (2003) as a lower limit. This 
value results from the (15 ± 4)% absorption detected over ~ 1 /3 
of the Hi Lyman a line. It is a lower limit because there is no 
way to measure the absorption over the central part of the line 
hidden by the interstellar medium absorption. Higher m values, 
(8^40) x 10 10 g s _1 , have recently been reported by Linsky et 
al. (2010). As discussed below, these authors have inferred hy- 
drogen mass-loss rates from a (8 ± 1)% absorption in the emis- 
sion line of singly ionized carbon (C 11) during the transit of the 
planet. 

Lecavelier des Etangs et al. (2010) have estimated the mass- 
loss rate of HD 189733b from a Hi Lya transit light curve, for 
different values of the X/EUV flux. For a best-fit value Fx/euv 
of20 times the solar value, they report a 1-cr range of 3.2xl0 9 < 
m< 1.4xlO n gs -1 . 

The X-ray luminosities are taken from Kashyap et al. (2008) 
for HD 209458 and Pillitteri et al. (2010) for HD 189733. A 
weakly efficient evaporation, with r] ~ 1%, seems to excluded 
from the observations of HD 189733. In contrast, only values 
of 77 > 1 are compatible with the observations of HD 209458b, 
which is physically not plausible. In fact, it is not clear what pos- 
sible energy source could significantly cumulate with the stellar 
radiation to drive such a "super-efficient" mass loss. The situa- 
tion is even more troublesome for HD 209458b if one considers 
the upper limit on Lx/euv given by Sanz-Forcada et al. (2010). 

The uncertainties on the measured values of Lx/euv and m 
are a first possible explanation to the apparent "super-efficiency" 
of HD 209458b's mass loss. In fact, a large part of the EUV 
spectrum (roughly from 200 to 912 A) does not contribute to 
the "X/EUV" flux measured with Rosat or XMM-Newton. If 
HD 209458 is as luminous between 200 and 912 A as it is 
below 200 A, then the position of the planet in the power di- 
agram (Fig. HJ would become compatible with a very efficient 
(77 ~ 100%) mass loss. However, in the case of the lower lumi- 
nosity estimation of Sanz-Forcada et al. (2010), the star would 
have to be orders of magnitude more luminous in the EUV above 
200 A than below. 

Another possible explanation for the super-efficiency of 
HD 209458b's mass loss may be found in the way the mass-loss 
rate was estimated by Linsky et al. (2010), which strongly de- 
pends on the assumed density profile of the exosphere, up to the 
Roche lobe (A. Lecavelier des Etangs, private communication). 

It is also possible that Lx/euv was larger during the observa- 
tions of Linsky et al. (2010) than in 2003. An increase by 1 or 2 
order(s) of magnitude, from the value of Kashyap et al. (2008) 
or Sanz-Forcada et al. (2010), respectively, is needed to have 
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Fig. 4. The power diagram for evaporating extrasolar planets. 
The "mass-loss power" L ln (Eq. [5j is plotted as a function of the 
X/EUV luminosity. Constant rj values are represented by the dot- 
ted lines. The hatched regions show the position of HD 189733b 
(orange) and HD 209458b (crimson). For HD 209458b, the 
right-pointing arrow show the upper limit on the X/EUV lumi- 
nosity reported by Sanz-Forcada et al. (2010). 



r\ values below 1, within the error bars. Such a large X lumi- 
nosity variation might not be expected for a solar-type star like 
HD 209458. Meanwhile, the only way to solve this question is 
to obtain simultaneous measurements of the X-ray flux and the 
escape rate. This should be possible with current space instru- 
mentation. 

Assuming the X/EUV luminosity from Kashyap et al. (2008) 
and the escape rate m between 10 10 and 4 x 10 11 g s _1 , it is pos- 
sible to infer exospheric properties out of the power diagram. 
In fact, while the m estimation of Vidal-Madjar et al. (2003), 
marginally compatible with r\ = 1, is based on the observation 
of an absorption in the H i stellar emission, the higher estimation 
of Linsky et al. (2010) is derived from the observation of an ab- 
sorption in the singly ionized carbon (C n) stellar emission. For 
this line, these authors measured an (8 + 1)% absorption, yielding 
m(C n) = 2. 1 x 10 7 g s . To go from this value to a total (hydro- 
gen) mass-loss rate m « 8 x 10 10 g s _1 , they rely on the assump- 
tions that (i) all the carbon is ionized (i.e., [C n/C] = 1) and (ii) 
that the carbon abundance is solar ([C/H] G = 2.7xl0~ 4 ; Asplund 
et al. 2009) in the upper atmosphere of HD 209458b. The total 
(hydrogen) mass-loss rate is then m = m(Hi) = m(Cn)/[C/H]. 

Therefore, an apparent super-efficient evaporation process 
with j] > 1, as shown in Fig. [4] could mean that [C/H] is un- 
derestimated. The [C/H] ratio would have to be increased by 
at least an order of magnitude with respect to [C /H] to give a 
m value compatible with 77 = 1, assuming Lx/euv values from 
Kashyap et al. (2008). On the other hand, if log Lx/euv 

< 26.12 

(Sanz-Forcada et al. 2010), then the upper atmosphere should be 
extremely carbon-rich with [C/H]/[C/H] G > 100. An underesti- 
mated [C/H] could also imply that the carbon ionization fraction 
is unity, as predicted by models of HD 209458b's atmosphere 
(Garcfa-Munoz 2007). 

Nevertheless, we caution that the measurements of m and 
Lx/euv for transiting planets may not be accurate enough at the 
present time to allow firm conclusions to be drawn out of the 
power diagram. Meanwhile, this shows how it could be possible 
to constrain the exospheric composition when precise measure- 
ments are available. 



5. Stability of evaporating atmospheres 

The lifetime of evaporating gaseous planets can be estimated as 
the time needed to exhaust the available reservoir of gas for a 
given mass-loss rate. An approximation of this life time can be 
obtained in an "integrated" version of the energy diagram shown 
in Fig.Q] We use the total potential energy E' p of a planet includ- 
ing tidal forces, defined by Eq. (10) in Lecavelier des Etangs 
(2007). 

The mean energy received (dLx/Euv/df) is calculated by in- 
tegrating Eq. (f2]i over a given time interval. Doing so, one should 
account for the evolution of Lx/euv(0 in time. In fact, the high- 
energy radiation output from a star is stronger when the star is 
younger and varies with time as a power-law Lx/Euv(f > T sat ) oc 
t~ a (e.g., Penz et al. 2008), with a close to 1. Stars younger 
than r sat are in the "saturation regime" where the evolution of 
Lx/euv(0 is almost flat. As discussed by Davis & Wheatley 
(2009), the duration r sat of the saturation period depends on the 
spectral type (see, e.g., Reiners et al. for M dwarfs). Here, we 
consider the equations used by Sanz-Forcada et al. (2010) 



^X/EUV 
^X/EUV 



= 6.3 x 10~ 4 L bol 
= 1.89 x 10 28 r L55 
= 2.03 x 10 20 Lbo°i 65 - 



(f < Tsat), 
(t > T sat ), 



(6) 

(7) 
(8) 



We use this set of equations to estimate (Lx/euv), the mean 
amount of X/EUV irradiation received during the life time of 
our sample stars per billion year, 



<L X /euv> = L X /Euv(f)d? + £ L X /Euv(f)dfj 

After injecting Eqs. (0 to © into Eq. (|9), we obtain 



(9) 



(Lx/euv) 



1[1.28 



x 10 17 L£ o 3 , 5 + 2.33 x 10 I7 L^ 575 



- 3.44 x 10 28 t 



where Lt, i and t+ are the present-day bolometric luminosity in 
erg s _1 and age of the star in Gyr, respectively. We calculate 
Lboi for our sample stars with Stefan-Boltzmann's law, Lboi = 
AncrTzJi 2 ^, where the effective temperatures and stellar radii are 
taken from the literature and listed in Table [TJ For M, K, G, F, 
and A stars, we assumed typical stellar ages of 10, 7, 5, 3, and 
1 Gyr, respectively. 

We want to know by how much the X/EUV irradiation is 
underestimated when the time evolution of Lx/euv is ignored. 
For this purpose, we introduce the factor y (see also Lecavelier 
des Etangs 2007), 



y 



(Lx/euv) 
Lx/euv 



(10) 



where Lx/euv at the denominator is calculated using Eq. (fTJ with 
t - t*. The values of y are listed in Table [TJ We found a mean 
value of y = 38 + 9 (lcr). We apply this y correction to the values 
of (d£x/Euv/df) calculated by integrating Eq. (0. 

This parameterization neglects the fact that the heating effi- 
ciency might be notably smaller at early times when the stellar 
luminosity is higher. Under high X/EUV illumination, Murray- 
Clay et al. (2009) have indeed suggested that the mass-loss 
regime should depart from an energy-limited escape (m oc 
^x/euv) anc ' i nstea d become "radiation/recombination"-limited 
(m oc F^myy). This effect would increase the lifetimes calcu- 
lated here, where we assume that the heating efficiency remains 
constant over time. 
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Fig. 5. The "integrated" energy diagram for 98 transiting planets: 
the energy needed to escape the whole planet mass versus the 
X/EUV flux received per billion year. The dotted lines represent 
constant lifetimes of 1, 5, and 10 Gyr (from bottom to top). Data 
points are calculated with rj = 0.15. Variations of r\ between 0.01 
and 1 are represented as horizontal gray error bars. 



These integrated energies are calculated for our sample of 
transiting planets and the result is shown in Fig. |5]for different 
values of tj. The ratio between -E' p (in erg) and (d£x/Euv/df) (in 
erg Gyr -1 ) is equivalent to a lifetime. As can be seen in Fig. [5] 
only a few planets have lifetimes below 10 Gyr. Meanwhile, this 
figure also highlights that any firm conclusions have to rely on 
an assumed value of rj. For instance, if 77 is close to 1 for WASP- 
17b, WASP-33b, or WASP-12b, then the atmospheres of these 
planets should be lost in less than 1 Gyr. Noteworthy, WASP- 
12b may well fill its Roche lobe fully enough that it is losing 
mass faster than the process considered here would predict (Li 
et al. 2010). 

An interesting case is that of Corot-7b. If this planet were 
made of gas, its whole mass would be evaporated within % 5 Gyr 
for 77 S 0.15. Considering that Corot-7b is probably made of rock 
(Leger et al. 2009, 201 1; Queloz et al. 2009), this means that the 
atmosphere of this planet could have been completely evapo- 
rated. This would make Corot-7b the first evaporation remnant 
detected. Possible origins of Corot-7b are discussed by Jackson 
et al. (2010). These authors conclude that this planet could have 
indeed started with a gas envelope, which total mass should have 
been < 2/3 M v 



6. Conclusion 

We have calculated mass-loss rates for transiting exoplanets us- 
ing an energy diagram. This approach assumes that mass loss 
occurs as a consequence of a massive stellar irradiation. Mass- 
loss rates for the complete sample of transiting planets are found 
between 10 6 and 10° g s _1 . This (large) range encompass all 
possible values of the atmospheric heating efficiency. 

We have shown that measurements of the heating efficiency 
can be constrained in a power diagram when the mass-loss rate, 
derived from observations of an evaporating atmospheres, and 
the X/EUV flux are both estimated. At the present time, it is 
only possible to spot two planets in this diagram. While the pre- 
cision on heating efficiency is limited by the estimation of the 
mass-loss rate, we have shown how it could be possible to infer 
exospheric properties, such as relative abundances, out of more 
precise measurements. In particular, simultaneous transit obser- 



vations in UV and X-rays, feasible with present-time space in- 
strumentation, will bring stronger constraints on the atmospheric 
heating efficiencies. 

The mass-loss rates are used to determine the lifetimes of 
evaporating planets. We have found that few planets have life- 
times below 1 Gyr but this requires a highly efficient heating 
with r] ~ 1 . 
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Table 1. Properties of known transiting planets and their host stars. 
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